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Abstract 

We report a measurement of the Rossiter-McLaughlin effect in the transiting extrasolar planetary sys- 
tem TrES-1, via simultaneous spectroscopic and photometric observations with the Subaru and MAGNUM 
telescopes. By modeling the radial velocity anomaly that was observed during a transit, we determine the 
sky-projected angle between the stellar spin axis and the planetary orbital axis to be A = 30 ±21 [deg]. 
This is the third case for which A has been measured in a transiting exoplanetary system, and the first 
demonstration that such measurements are possible for relatively faint host stars {V ^ 12, as compared to 
y ^ 8 for the other systems). We also derive a time of mid-transit, constraints on the eccentricity of the 
TrES-lb orbit (e = 0.048 ± 0.025), and upper limits on the mass of the Trojan companions (^14 M®) at 
the 3ct level. 

Key words: stars: planetary systems: individual (TrES-1) — stars: rotation — techniques: photo- 
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1. Introduction 

The Rossiter-McLaughlin effect (hereafter the RM ef- 
fect) is a phenomenon originally reported as a "rotational 
effect" in eclipsing binary systems by Rossiter (1924) (for 
the Beta Lyrae system) and McLaughlin (1924) (for the 
Algol system) . In the context of extrasolar planetary sci- 
ence, the RM effect is seen as a radial velocity anomaly 
during a planetary transit caused by the partial occul- 
tation of the rotating stellar disk (see Ohta et al. 2005, 
Gimenez 2006, or Gaudi & Winn 2007, for theoretical de- 
scriptions). The radial velocity anomaly depends on the 
trajectory of the planet across the disk of the host star, 
and in particular on the alignment between that trajec- 
tory and the rotation field of the star. By monitoring this 
anomaly throughout a transit one can determine whether 
or not the planetary orbital axis is well-aligned with the 
stellar spin axis. In the Solar System, the orbits of all 8 
planets are known to be well-aligned with the solar equa- 

* Based in part on data collected at Subaru Telescope, which is 

operated by the National Astronomical Observatory of Japan. 
** JSPS Fellow. 



tor, but this is not necessarily the case for exoplanetary 
systems, or for hot Jupiters in particular. The key pa- 
rameter is the sky-projected angle between the stellar spin 
axis and the planetary orbital axis. A, and measurements 
of this "misalignment angle" for various exoplanetary sys- 
tems will help to place the Solar System in a broader con- 
text. 

Specifically, measurements of the RM effect for exoplan- 
etary systems are important because of the implications 
for theories of migration and hot Jupiter formation. So 
far, measurements of A for two systems have been re- 
ported; Queloz et al. (2000) and Winn et al. (2005) for 
HD 209458, Winn et al. (2006c) for HD 189733. In both 
of those systems, the host star is very bright (^ ~ 8), 
facilitating the measurement. The observed values of A 
for the two systems are small or consistent with zero, 
which would imply that the standard migration mech- 
anism (planet-disk interaction) does not alter the spin- 
orbit alignment grossly during the planetary formation 
epoch. However, just these few examples are not enough 
for statistical constraints on other hot Jupiter forma- 
tion theories, including planet-planet interaction (Rasio & 
Ford 1996; Weidenschilling & Marzari 1996), the "jumping 
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Jupiter" model (Marzari & Weidenschilling 2002), or the 
Kozai mechanism (Wu & Murray 2003), which may lead 
hot Jupiters to have significantly misaligned orbits. Thus 
further measurements of the RM effect for other transiting 
systems are valuable. Given that most of ongoing tran- 
sit surveys target relatively faint {V ~ 12) host stars, it 
is also important to extend the reach of this technique to 
fainter stars. Further observations for new targets would 
be useful to constrain planet formation theories, and more 
importantly, have a potential to discover large spin-orbit 
misalignments, which would be a challenge to some theo- 
retical models. 

In this paper, we report the measurement of the RM ef- 
fect and the constraint on A for TrES-1 (V = 11.8) which 
has a significantly fainter host star than those in previous 
studies (y ~ 8). In addition to the fainter host star, this 
work differs from previous studies of the RM effect in that 
we have conducted simultaneous spectroscopic and photo- 
metric observations. This new strategy offers several po- 
tentially important advantages. First, the simultaneous 
photometry eliminates any uncertainty in the results due 
to the orbital ephemeris and the transit depth. Although 
this did not turn out to be crucial for the present work, 
it will be useful for newly discovered targets which still 
have uncertainty in the times of transits. Second, the 
transit depth might be expected to vary due to star spots 
or transient events, and indeed evidence for star spots 
was reported in HST/ACS photometry for this system 
(Charbonneau et al. 2007). Thus simultaneous monitoring 
is useful to assess anomalies in the transit depth. Finally, 
obtaining all of the data on a single night is useful to 
avoid systematic errors in radial velocities from long-term 
instrumental instabilities. Moreover, although it need not 
be simultaneous, the photometry also helps to determine 
the limb-darkening parameter for the visual band, which 
can be used in the interpretation of the RM-affected spec- 
tra. In this way, we can determine the limb-darkening 
parameter of the host star directly from the data, instead 
of assuming a value based on stellar atmosphere models. 

We describe our observations in Section 2 and report 
the results in Section 3. Section 4 provides a discussion 
and summary of this paper. 

2. Observations and Data Reduction 

We observed the planet-hosting KOV star TrES-1 on UT 
2006 June 21, the night of a predicted planetary transit, 
using the Subaru 8.2 m Telescope at Mauna Kea and the 
MAGNUM 2 m Telescope at Haleakala, both in Hawaii. 
The event is predicted as the 238th transit from the first 
discovery, namely E = 238 {E: integer) in the ephemeris 
by Alonso et al. (2004); 

Tc{E) = 2453186.8060(±0.0002) 

+ Ex 3.030065(±0.000008) [HJD]. (1) 

The transit occurred shortly after midnight. We observed 
TrES-1 during 5 hours bracketing the predicted transit 
time, through air masses ranging from 1.0 to 1.3. 



2.1. MA GNUM Photometry 

The MAGNUM 2 m telescope is located near the 
Haleakala summit on the Hawaiian Island of Maui 
(Kobayashi et al. 1998; Yoshu 2002; Yoshu et al. 2003). 
The MAGNUM photometric observation was conducted 
in parallel with the Subaru spectroscopic observation. We 
employed the Multi-color Imaging Photometer (MIP) us- 
ing a 1024 X 1024 pixel CCD with a V band filter, covering 
4750 A < A < 6180 A, and we set 9 dithering positions 
(3x3 positions) on the CCD. The MIP has a 1.5 x 1.5 
arcmin^ field of view (FOV) with a pixel scale of 0.277 
arcsec/pixel. We used 2MASS J19041058-H3638409 as our 
comparison star for differential photometry. This star is 
close enough to fit in the MIP field of view, and is known 
to be photometrically stable at a level sufficient for our 
study (e.g., Charbonneau et al. 2005). The exposure time 
was either 40 or 60 seconds according to observing condi- 
tions so that the photon counts are close to the saturation 
level of the CCD, with a readout/setup time of 60 seconds. 

We then reduced the images with the standard MIP 
pipeline described in Minezaki et al. (2004). We deter- 
mined the apparent magnitudes of TrES-1 and the com- 
parison star using an aperture radius of 20 pixels. The 
typical FWHM of each star ranged from 1.4 to 1.9 arcsec 
(from 5 to 7 pixels). We estimated the sky background 
level with an annulus from 20 to 25 pixels in radius cen- 
tered on each star, and subtracted the estimated sky con- 
tribution from the aperture flux. Then we computed the 
differential magnitude between TrES-1 and the compari- 
son star. After these steps, we decorrelated the differential 
magnitude from the dithering positions and eliminated ap- 
parent outliers from the light curve, most of which were 
obtained at the 9th dithering position. We do not find 
any clear correlations with other observing parameters. 

For the analysis of transit photometry, Pont et al. 
(2006) studied the time-correlated noise (the so-called 
"red noise") in detail, and Gillon et al. (2006) introduced 
a simple and useful method to account at least approx- 
imately for the effect of the red noise on parameter es- 
timation. Based on these studies, we used the following 
procedure to determine the appropriate data weights for 
the MAGNUM photometry (which are similar to that em- 
ployed by Winn et al. 2006a) . 

We first fitted the MAGNUM light curve with the ana- 
lytic formula given in Ohta et al. (2005) and foimd the 
residuals between the data and the best-fitting model. 
Using only the Poisson noise as an estimate of the er- 
ror in each photometric sample, we found x^/^dof ~ 2.8 
[vdof'- degrees of freedom), implying that the true errors 
are significantly in excess of the Poisson noise. We also 
found the residuals in the early part of the night (before 
^ 2453907.91 [HJD]) to be significantly larger than those 
from later in the night. This larger scatter could have 
been caused by shaking of the telescope by the stronger 
winds that occurred during the early part of the night. 
We thus estimated error-bars separately for the early part 
of the night and the later part of the night, as described 
below. 
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Fig. 1. Top: A photometric light curve from the MAGNUM 
observation (184 samples). The error-bars are scaled to satisfy 
X^/^dof = 1.0 (see Section 2.1.). Bottom: 20 radial velocity 
samples computed from the Subaru/HDS spectra. The values 
and uncertainties are presented in Table 1. 



Table 1. Radial velocities obtained with the Subaru/HDS. 



Time [HJD] Value [m s '-] Error [m s 



2453907.87018 


18.7 


14.0 


2453907.88139 


30.5 


12.5 


2453907.89262 


54.6" 


12.0 


2453907.90384 


24.3 


10.4 


2453907.91506 


26.4 


11.4 


2453907.92628 


30.4 


10.9 


2453907.93750 


22.4 


14.3 


2453907.94873 


2.9 


11.0 


2453907.95996 


-7.1 


12.1 


2453907.97119 


-22.3 


13.3 


2453907.98241 


-40.5 


13.3 


2453907.99364 


-39.2 


13.0 


2453908.00488 


-9.8 


12.2 


2453908.01610 


-30.5 


13.8 


2453908.02732 


-17.7 


13.6 


2453908.03854 


-24.7 


12.2 


2453908.04978 


-27.5 


11.1 


2453908.06100 


-38.2 


13.3 


2453908.07223 


-23.7 


11.2 


2453908.08345 


-23.0 


9.6 



^: A possible outlier. 



First, we rescaled the error bars to satisfy 'yd Ivdof = 1-0 
(step 1), namely 0.00259 for the early data and 0.00189 
for the late data. The light curve with these rescaled 
error bars is shown in the upper panel of Fig. 1. Next, 
in order to assess the size of time-correlated noise for the 
MAGNUM data, we solved the following equations, 

o\=oI + gI, (2) 

where (Ti is the standard deviation of each residual and 
is the standard deviation of the average of the successive 
N points, aw is called the white noise, which is uncor- 
related noise that averages down as while a,, is 

called the red noise, which represents correlated noise that 
remains constant for specified N. We calculated a^j and 
(7r for the choice A'^ = 30 (corresponding to one hour), find- 
ing a.u] = 0.00219 and ar = 0.00139 for the early data. On 
the other hand, we found cr2 < (suggesting a smaller level 
of the red noise) for the late data. The choice of vV = 30 
is fairly arbitrary; other choices of N between 5 and 50 
gave similar results. We then adjusted the error bars for 
the early night by multiplying [1 -|- iV {(Tr/cFw)^]^^^ 3.6 
(step 2). We did not change the error bars for the late 
night data. We adopted these rescaled uncertainties for 
subsequent fitting procedures. 

2.2. Subaru HDS Spectroscopy 

We used the High Dispersion Spectrograph (HDS) on 
the Subaru Telescope (Noguchi et al. 2002). We employed 
the standard I2a set-up of the HDS, covering 4940 A < 
A < 6180 A with the Iodine absorption cell for measuring 
radial velocities. The slit width of 0"8 yielded a spectral 



resolution of ~ 45000, and the seeing was between 0''75 
and 1''2. The exposure time for TrES-1 was 15 minutes 
yielding a typical signal-to-noise ratio (SNR) ~ 60 per 
pixel. In order to estimate systematic errors from short 
term instrumental variations, we also obtained spectra of 
the much brighter {V = 4.7) KOV star HD 185144 before 
and after the series of TrES-1 exposures. This star is 
known to be stable in velocity (Johnson et al. 2006). We 
obtained five 30 s exposures of HD 185144, each having a 
SNR of approximately 200 pixel^^. 

We processed the frames with standard IRAF^ proce- 
dures and extracted one-dimensional spectra. Next, we 
calculated relative radial velocity variations by the algo- 
rithm following Sato et al. (2002). We used this algo- 
rithm because it properly takes into account the fairly 
large changes of the instrumental profile during the obser- 
vations. The HDS is known to experience appreciable in- 
strumental variations even within a single night, reported 
in Winn et al. (2004) and Narita et al. (2005). We es- 
timated internal errors of the radial velocities from the 
scatter of the radial velocity solutions for 2A segments 
of the spectra. The typical errors are 10 ~ 15 [m s^^], 
which are reasonable values to be expected from the pho- 
ton noise limit. Note that we do not find any evidence 
of star spots or transient events during our photometric 
observation (see Fig. 1). For this reason, we have not 
accounted for possible systematic errors in the velocities 
due to star spots. We also reduced the HD 185144 spectra 

^ The Image Reduction and Analysis Facility (IRAF) is dis- 
tributed by the U.S. National Optical Astronomy Observatories, 
which are operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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with the same method in order to check for systematic er- 
rors due to short-term instrumental instabilities. The rms 
of the radial velocity of HD 185144 is less than 5 [m s~^], 
attesting to good instrumental stability. The resultant 
radial velocities of TrES-1 and their errors are shown in 
Table 1 and the lower panel of Fig. 1. 

3. Results 

As described above, we have obtained 20 radial veloc- 
ity samples and 184 V band photometric samples taken 
simultaneously covering the transit. In addition, in or- 
der to search for an optimal solution of orbital param- 
eters for TrES-1, we incorporate our new data with 12 
previously published radial velocity measurements using 
the Keck I telescope (7 by Alonso et al. 2004 and 5 by 
Laughlin et al. 2005) and 1149 z band photometric mea- 
surements spanning 3 transits using the FLWO 1.2m tele- 
scope (Winn et al. 2006b). The uncertainties of the FLWO 
data had already been rescaled by the authors such that 
X^/'^dof = 1.0 for each transit (namely, the step 1 has 
been done). For the step 2, we find < for these 
data, thus we did not modify these error bars further. 
We employ the analytic formulas of radial velocity and 
photometry including the RM effect given in Ohta et al. 
(2005) and Ohta ct al. (2006) (hereafter the OTS formu- 
lae) in order to model the observed data. Note that based 
on the previous studies by Winn et al. (2005) and Winn 
et al. (2006c), we have learned that the OTS formulae 
systematically underestimate the amplitude of radial ve- 
locity anomaly by approximately 10%. This is possibly 
because the radial velocity anomaly defined by Ohta et al. 
(2005) and that measured by the analysis pipeline are dif- 
ferent. Thus we correct V sin Is in the OTS formulas by 
modifying ysin/s(OTS) = F sin (real) * 1.1. This cor- 
rection presumably gives more realistic values for V sin Is 
and A, and has little influence on any of the other pa- 
rameters. Here we assume circular orbits of the star and 
the planet about the center of mass (namely, e = 0) . We 
adopt the stellar mass Ms = 0.87(±0.03) [Mq] (Laughlin 
et al. 2005) and the orbital period P = 3.0300737 [days] 
and Tc(0) = 2453186.80603 [HJD] (Winn et al. 2006b). 
As a result, our model has 15 free parameters in total. 
Eight parameters for the TrES-1 system include the ra- 
dial velocity amplitude K, the sky-projected stellar rota- 
tional velocity V sin Is, the misalignment angle between 
the stellar spin and the planetary orbit axes A, the lin- 
ear limb-darkening parameter for V band uy, the same 
for z band Uz, the ratio of star-planet radii Rp/Rg, the 
stellar radius Rs, and the orbital inclination i. Here we 
assume that the limb-darkening parameters for the spec- 
troscopic and photometric models are the same.^ We also 



^ This is likely to be a good approximation because the photomet- 
ric band is a good match to the region with abundant I2 lines 
where the radial velocities are measured. However, in princi- 
ple the correspondence is not exact because the limb-darkening 
function may not be identical in the lines as opposed to the 
continuum, and because the influence of limb-darkening on the 
RV-measuring algorithm has yet to be investigated in detail. 



add three parameters for velocity offsets to the respective 
radial velocity dataset vi (for our template spectrum), 
V2 (for Alonso et al. 2004) and V3 (for Laughlin et al. 
2005), and four parameters for the times of mid-transit 
Tc{E) [E = 234,235,236,238). 

In previous studies of the exoplanetary RM effect, it 
was possible and desirable to determine both Vsinlg and 
A from the radial velocity data. In this case, there are 
two reasons to prefer an external determination of l^sin/g. 
First, the signal-to-noisc ratio of the anomaly is smaller, 
because of the faintness of the host star. This makes it 
valuable to reduce the number of degrees of freedom in 
the model. Second, the transit geometry is nearly equato- 
rial, which introduces a very strong degeneracy between 
VsinJs and A, as explained by Gaudi & Winn (2007). The 
alternative we have chosen is to adopt a value for V sin Is 
based on previous observations, and use the radial- velocity 
anomaly to determine A. (We have also investigated our 
ability to determine both parameters, as described be- 
low.) Laughlin et al. (2005) reported Vsinis = 1.08±0.30 
[km s^^] for the TrES-1 host star from their analysis of 
the observed spectral line profiles; this is the most reliable 
estimate for Vsinis to date. We incorporate this infor- 



■ Vsin/s-1.081 



mation into our model by adding a term |^ ^ 
the fitting statistic. Thus our statistic is 



to 



E 



'^i.obs '^i.calc 



2 Ary = 1333 

+ E 



fj,obs /- 



j^calc 



Vsinis - 1.08 



0.30 



where Vcaic and fcaic represent the values calculated by the 
OTS formulae with the above parameters. We find opti- 
mal parameters by minimizing the statistic of Eq. (4) 
using the AMOEBA algorithm (Press et al. 1992), and es- 
timate confidence levels of the parameters using Ax^ from 
the optimal parameter set. To assess the dependence of 
our results on the a priori constraint on Vsinlg, wc also 
compute and compare the best-fit values and uncertainties 
by using another function: 



X 



E 



Ar^ = 1333 

- E 



j^obs fj.calc 



for reference. In addition, we note that the third radial 
velocity sample of our data (t = 2453907.89262 [HJD] ) may 
appear to be an outlier, but it lies just about Sa from a 
theoretical radial velocity curve (e.g.. Fig. 2). For clarity, 
we calculate Eq. (4) and Eq. (5) with and without that 
sample. 

The results for both x^ statistics are presented in 
Table 2. The minimum is 1308.57 (1296.39) for Eq. 
(4) and 1305.18 (1298.57) for Eq. (5) with 1350 (1349) 
degrees of freedom, where the numbers in parentheses re- 
fer to the case without the outlier. In Fig. 2, we present 
the radial velocities and the best-fit curve with (the top 
figure, marked with "a") and without (the middle figure, 
marked with "b") the a priori constraint. In addition. 
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Fig. 3. Contours of constant in (Vsin/s, A) space, based on simultaneous fitting of 32 radial velocity samples and 1333 photometric 
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Ax^ = 2.30 (inner) and Ax^ = 6.17 (outer), while the dotted line shows Ax^ = 1.00 (inner) and Ax^ = 4.00 (outer), respectively. 
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Table 2. Best-fit values and uncertainties^ of the free parameters. 





AU data 
With the constraint 


All data 
Without the constraint 


Subaru/MAGNUM 
With the constraint 


Parameter 


Value 


Uncertainty 


Value 


Uncertainty 


Value 


Uncertainty 


(All rv samples) 
K [m s-i] 


113.1 


±2.5 


113.1 


±2.5 


115.21= 


fixed 


Vsvals [km s^-"-] 


1.3 


±0.3 


2.5 


±0.8 


1.3 


±0.3 


A [deg] 


30 


±21 


48 


±17 


28 


±24 


My 


0.57 


±0.05 


0.57 


±0.05 


0.59 


±0.05 




0.37 


±0.03 


0.37 


±0.03 


— 


_ 


Rp/Rs 


0.1382 


±0.006 


0.1382 


±0.006 


0.13686^= 


fixed 


Rs [Re] 


0.82 


±0.02 


0.82 


±0.02 


0.8ir 


fixed 


i [deg] 


88.4 


±0.3 


88.4 


±0.4 


88.9^= 


fixed 


vi [m s""'"] 


1.3 


±3.0 


4.0 


±3.5 


0.7 


±2.9 


V2 [m s"i] 


-0.2 


±4.9 


-0.2 


±4.9 


- 


— 


W3 [m s"^] 


-5.5 


±1.6 


-5.5 


±1.6 


- 


- 


rc(234)- 2453000 [HJD] 


895.84298 


±0.00015 


895.84298 


±0.00015 






Tc(235) - 2453000 [HJD] 


898.87342 


±0.00014 


898.87342 


±0.00014 






Tc(236) - 2453000 [HJD] 


901.90371 


±0.00016 


901.90371 


±0.00016 






Tc(238) - 2453000 [HJD] 


907.96407 


±0.00034 


907.96407 


±0.00034 


907.96408 


±0.00034 


(Without the outlier) 
K [m s-i] 


112.8 


±2.5 


112.8 


±2.5 


115.2'= 


fixed 


Vs'mis [km s^-'] 


1.3 


±0.3 


2.1 


±0.8 


1.3 


±0.3 


A [deg] 


24 


±23 


39 


±21 


22 


±27 


vi [m s ^] 


1.1 


±3.1 


0.9 


±3.6 


-1.5 


±3.0 



"iComputed by Ax^ = 1.00. ''lAlonso et al. (2004). ^^iWinn et al. (2006b). 



we also compute the best- fit curve without the constraint, 
but assuming that A = [deg] (the bottom figure, marked 
with "A = [deg]"). Fig. 3 plots {VsinIs,X) contours cal- 
culated with (the left panel, marked with "a" ) and without 
(the right panel, marked with "b") the a priori constraint. 
Note that we only show here the results with the possi- 
ble outlier in Fig. 2 and Fig. 3, since the same figures 
but without the outlier have basically similar appearance 
and have less information to show. As a result, we find 
A = 30±21 (24±23) [deg] and V^sin/, = 1.3±0.3 (1.3±0.3) 
[km s^'^] for TrES-1 with the a priori constraint, and our 
findings except for A are in good agreement with previous 
studies (Alonso et al. 2004; Laughhn et al. 2005; Winn 
et al. 2006b). On the other hand, the result without 
the a priori constraint, A = 48 ±17 (39 ±21) [deg] and 
Fsin/s = 2.5 ±0.8 (2.1 ±0.8) [km s^^], agrees with the 
above result within about la for A and Vsinlg, and is 
also consistent with the previous results for other param- 
eters. In case we calculate without the constraint but 
assuming that A = [deg] , we find that the minimum 
is 1309.85 (1298.91) with 1351 (1350) degrees of freedom, 
and Vs'mis ~ 1.5 ±0.5 (1.5 ±0.5) [km s^^]. Consequently, 
our results for A have fairly large uncertainties. We find 
at least that the orbital motion of TrES-lb is prograde. 
Additional radial velocity measurements during transits, 
and a more precise measurement of l^sin/g, would be de- 
sirable to pin down A and help to discriminate between 
different modes of migration. 



4. Discussion and Summary 

We have presented simultaneous spectroscopy and pho- 
tometry of a transit of TrES-lb, exhibiting a clear de- 
tection of the Rossiter-McLaughlin effect and consequent 
constraints on the alignment angle between the stellar spin 
and the planetary orbital axes. Our philosophy has been 
to use all of the best data available at present. However, it 
is also interesting to examine how well we are able to deter- 
mine the system parameters using only the data gathered 
on a single night using Subaru and MAGNUM. This is be- 
cause for future studies of newly discovered transiting ex- 
oplanetary systems, higher-precision data from other ob- 
servatories may not be available. We repeat the fitting 
procedure without the Keck and FLWO data but still as- 
suming the system parameters other than Tasinis, A, uy, 
vi, and Tc(238) to be the values presented in Alonso et al. 
(2004) and Winn et al. (2006b). We find almost the same 
values and uncertainties for the parameters above (the 
right side of Table 2) as before, indicating that a single 
night's data would have done almost as well as the full 
data set. 

Using transit timing of the TrES-1 system, Steffen & 
Agol (2005) reported a constraint on the existence of ad- 
ditional planets in the system. Subsequently, Winn et al. 
(2006b) pointed out that Tc{E) {E = 234,235,236) are con- 
sistent with their ephemeris at about the 2a level, but oc- 
curred progressively later than expected. Here we present 
Tc{E) {E = 238) in Fig. 4. Our result is also consistent 
with the published ephemeris within 1.5cr, occurring only 
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Fig. 4. Timing residuals of the time of mid-transit for 
each epoch (E), based on the ephemeris of Winn et al. 
(2006b). The results for E = 234,235,236 are determined 
by the FLWO photometry, while E = 238 is computed by 
the MAGNUM data. 



Fig. 5. Contours of constant i^i (ecosii;,esina;) space, 
based on simultaneous fitting of 32 radial velocity samples 
and 1333 photometric samples. At each grid point, all of 
the parameters except for {e,u!) were optimized. Contour 
line types are the same as those in Fig. 3 



slightly later than expected from Winn et al. (2006b). 

In addition to producing transit timing variations, any 
additional bodies in the TrES-1 system could excite the 
orbital eccentricity of the TrES-lb planet, which could be 
detectable in the radial velocity measurements of the host 
star. Thus it is worthwhile to put empirical constraints on 
the orbital eccentricity. For this purpose, we compute 
using Eq. (4) for fixed values of {e,uj) over numbers of grid 
points to map out the allowed region in (e, uj) space. In 
fact the most appropriate parameters are ecosw and esino; 
since their uncertainties are uncorrelated (see Fig. 5). The 
resulting constraints are ecosa; = 0.005±0.005 and esinw = 
—0.048 ± 0.024. Our result for ecosw is similar to that of 
Charbonncau et al. (2005), which was based on the timing 
of the secondary eclipse. The value of ecosw is consistent 
with zero within Icr, and the value of esinw is consistent 
with zero within 2a. Thus we do not find any strong 
evidence for a nonzero orbital eccentricity in the TrES-1 
system. 

There is another interesting application of our data. 
Recently, Ford & Gaudi (2006) studied the detectabil- 
ity of "hot Trojan" companions near the L4/L5 points 
of transiting hot Jupiters, through any observed differ- 
ence between the time of vanishing stellar radial veloc- 
ity variation (Tq) SiJid the time of the midpoint of the 
photometric transit (Tc). Our strategy of simultaneous 
spectroscopic and photometric transit observations is ide- 
ally suited for searching for the hot Trojan companions. 
For the TrES-1 system. Ford & Gaudi (2006) set an up- 
per limit on the mass of the Trojan companions ~ 51Mq 
at the 3f7 level (assuming the circular orbit), using the 
radial velocity samples of Alonso et al. (2004). Here 
we compute At = To — rc(238) using all available out-of- 
transit radial velocity samples with (without) the possi- 
ble outlier, and both for the circular and the eccentric 
orbit. We find At = 3.2 ± 11.2 (3.2 ± 11.8) [min] (circu- 
lar) and At = 33.0 ±52.6 (33.7 ±44.8) [min] (eccentric). 
Accordingly, we set constraints on the mass of the Trojan 



companions, Mt, which is defined as the difference in the 
mass at L4 (Mt,l4) and the mass at L5 (Mt.ls) (namely, 
AIt = Mt,l4 — Mt,l5), through the relation; 

47r At 

Mr^^M,^^, (6) 

where ctmt and cta* indicate the uncertainties of Mt and 
At, respectively (Eq. (1) and Eq. (2) in Ford & Gaudi 
2006). Note that we adopt Mp = 0.73 ±0.03 [Mjup] which 
is determined by this work. We find Mr = 1.2 ±4.3 (1.2 ± 
4.6) [Me] (circular) and Mt = 13 ±20 (13 ±17) [M®] (ec- 
centric). As a result, we exclude the Trojan companions 
near the L4 point more massive than ~ 14 (15) [M®] if the 
orbit is circular, and ~ 74 (65) [M^] if the orbit is allowed 
to be eccentric, both at the 3a level. Our constraint un- 
der the reasonable assumption of a circular orbit is more 
stringent than that of Ford & Gaudi (2006) by a factor of 
4, because we have increased the number of radial velocity 
samples and because our data cover the critical phase to 
determine Tq. Consequently, we conclude that we do not 
find any sign of the existence of additional bodies in the 
TrES-1 system at present. 

In this paper, we have placed a constraint on the sky- 
projected angle between the stellar spin axis and the 
planetary orbital axis for the TrES-1 system, namely 
A = 30 ±21 [deg] using all available data and informa- 
tion from previous studies. Although we can not dis- 
criminate whether the spin-orbit angle in this system is 
well- aligned or not at this point, our constraint on A 
clearly indicates the prograde orbital motion of TrES-lb. 
The uncertainty is larger than in previous studies (~ 1 
[deg] for HD 209458 and HD 189733) because the host 
star is significantly fainter in this case. Although fur- 
ther radial velocity measurements during transit would 
be necessary to pin down A more stringently, we have 
demonstrated for the first time that such measurements 
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are possible for such a faint target. This is important 
because most of the newly discovered transiting planets 
from ongoing transit surveys will have relatively faint host 
stars. For example, the new targets that were discovered 
in 2006, namely XO-1 (McCuUough et al. 2006), TrES-2 
(O'Donovan et al. 2006), HAT-P-1 (Bakos et al. 2006), 
WASP-1 and WASP-2 (Cameron et al. 2007), are all in 
this category. Combining future measurements of A in 
other transiting systems, we would be able to determine 
the distribution of A for exoplanetary systems with useful 
statistical accuracy. 
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